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Abstract 

A measurement of the underlying event (UE) activity in proton-proton collisions is 
performed using events with charged-particle jets produced in the central pseudora¬ 
pidity region (\rf et \ < 2) and with transverse momentum 1 < Pj 1 < 100GeV. The 
analysis uses a data sample collected at a centre-of-mass energy of 2.76 TeV with the 
CMS experiment at the LHC. The UE activity is measured as a function of Pj 1 in terms 
of the average multiplicity and scalar sum of transverse momenta (pj) of charged par¬ 
ticles, with \rj | < 2 and pj > 0.5 GeV, in the azimuthal region transverse to the highest 
pj jet direction. By further dividing the transverse region into two regions of smaller 
and larger activity, various components of the UE activity are separated. The mea¬ 
surements are compared to previous results at 0.9 and 7 TeV, and to predictions of 
several Monte Carlo event generators, providing constraints on the modelling of the 
UE dynamics. 
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1 Introduction 

Hadron production in high-energy proton-proton (pp) collisions originates from multiple scat¬ 
terings of the partonic constituents of the protons at central rapidities, and from "spectator" 
(noncolliding) partons emitted in the very forward direction. The produced partons reduce 
their virtuality through gluon radiation and quark-antiquark splittings, and finally fragment 
into hadrons at scales approaching 0.2 GeV (Aqcd)- Usually, one separates the produced hadrons 
into two classes: those coming directly from the fragmentation of partons resulting from the 
scattering with the largest momentum transfer (hard scattering) in the event, and the rest (un¬ 
derlying event, or UE). The UE thus consists of hadrons coming from (i) initial- and final-state 
radiation (ISR, FSR) from the hard scattering, (ii) softer partonic scatters in the same pp collision 
(multiple parton interactions, or MPI) possibly with their own initial- and final-state radiation, 
and (iii) proton remnants concentrated along the beam direction. 

An accurate understanding of the UE is required for precise measurements of standard model 
processes at high energies and searches for new physics. Indeed, the UE affects measure¬ 
ments of isolated high transverse momentum pj leptons or photons, and it dominates most 
of the hadronic activity from the overlapping pp collisions taking place in a given bunch cross¬ 
ing (pileup) at the high luminosities achieved by the CERN LHC. The semi-hard and low- 
momentum partonic processes, which dominate the UE, cannot be adequately calculated with 
perturbative Quantum Chromodynamics (pQCD) methods alone, and require a phenomeno¬ 
logical description containing parameters that must be tuned to data. 

The topological structure of pp interactions with a hard scattering can be used to define ex¬ 
perimental observables sensitive to the UE. One example is the study of particle properties in 
regions away from the direction of the products of the hard scattering. At the Tevatron, the 
CDF experiment measured UE observables using inclusive jet and Drell-Yan (DY) events in 
pp collisions at centre-of-mass energies Qs = 0.63, 1.8, and 1.96TeV [[l]-[3). In pp collisions 
at the LHC, the ALICE, ATLAS, and CMS experiments have carried out UE measurements at 
y/s = 0.9 and 7 TeV using events containing a leading (highest pj) charged-particle jet EMU or 
a leading charged particle EHE), or a DY lepton pair 0. In this paper, we study the UE activity 
in pp collisions at y/s = 2.76 TeV by measuring the average multiplicity and scalar transverse 
momentum sum (Zpj) densities of charged particles in the azimuthal region orthogonal to the 
direction of the leading charged-particle jet, referred to as the transverse region. 

At a given centre-of-mass energy, the UE activity is expected to increase with the momentum 
transfer between the interacting partons (hard scale). On average, increasingly hard parton 
interactions result from pp collisions with decreasing impact parameters between the two pro¬ 
tons, which in turn enhance the overall hadronic activity originating from MPI until a satu¬ 
ration is reached for central collisions with maximum overlap IfTOlUT l. At the same time, the 
activity related to the ISR and FSR components also increases with the hard scale. For events 
with the same hard scale, probed by the pj of jets or DY pairs, the MPI activity rises with yj s, 
as more partons are expected in the protons at increasingly smaller parton fractional momenta 
x ~ 2 pj/ y/s CUB EH. Hence, studying the UE as a function of the hard scale at several centre- 
of-mass energies provides an insight into the UE dynamics and its evolution with the collision 
energy, further constraining the model parameters. 

The paper is organised as follows. Section 2 presents the main features of the Monte Carlo 
(MC) event generators used in this study to provide a description of the UE properties. Section 
3 briefly describes the experimental methods, observables, event and track selection, as well as 
the corrections and systematic uncertainties of the measurements. The results are presented in 
Section 4, and summarised in Section 5. 
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2 Monte Carlo event generators 

In this analysis, the PYTHIA6 Ifl2l , PYTHIA8 [|T3| , and HERWIG++ Cl4l MC event generators are 
used with various tunes that are described below. In PYTHIA, the 2 —>2 parton scatterings, in¬ 
cluding MPI, are described by leading-order pQCD, with the 1 / p.\ cross section divergence reg¬ 
ularised by introducing a low-px infrared cutoff (pio), such that the diverging term is replaced 
by 1/ (pj + pj 0 ) 2 . There are various tunable parameters that control the behaviour of this reg- 
ularisation, the matter distribution of partons in the transverse plane within the hadrons, and 
the final-state colour reconnection effects among the produced partons. When QCD radiation 
is modelled via a px~ordered evolution, the MPI and parton showers are interleaved in one 
common sequence of decreasing pj values Itl5| . For the latest version of PYTHIA6 only ISR 
showers and MPI are interleaved, while in PYTHIA8 FSR showers are also included. The final 
nonperturbative transition of partons to hadrons is described by the Lund string fragmentation 
model Ifl6l . 

Another general-purpose generator, HERWIG++, is similar to PYTHIA, but uses angular-ordered 
parton showers and the cluster model B4l for hadronisation. It has an MPI model similar to 
the one used by PYTHIA, with tunable parameters for regularising the partonic cross sections 
at low momentum transfer, but does not include the interleaved evolution with ISR and FSR. 

Both MC models incorporate multiple parton collisions "perturbatively"—i.e. based on a "reg- 
ularisation" of the underlying pQCD subprocesses' cross sections — but require a nonpertur¬ 
bative ansatz for the impact parameter profile of the colliding protons. The frequency of MPI is 
then generated by assuming a Poissonian distribution of the number of elementary partonic in¬ 
teractions over the overlapping pp volume, with the average number depending on the impact 
parameter of the hadronic collision IllOlITTII . The MPI cross section is dominated by scatterings 
with semi-hard momentum transfers, 0(1-2 GeV), involving low-x partons, and thus shows a 
stronger dependence on the evolution of the low-px infrared cutoff, and on the incoming par- 
ton densities than the single hard-scattering interactions IflOl fill. In PYTHIA6, PYTHIA8 and 
HERWIG++, the energy dependence of MPI is mostly controlled by the energy evolution of 
the low-px infrared cutoff parameter ,which follows a (tunable) power law dependence on the 
centre-of-mass energy fl2til4ll . The UE activity accompanying various types of hard scattering 
processes is well described by MC event generators, ® [5, [ZHSl, illustrating the universality of 
MPI in different event topologies and hard-scattering production processes. Such a universality 
is confirmed by the similarity between the UE activity measured in DY f9l and jet-dominated 
events SIS, 7113, despite their different underlying parton radiation patterns. 

In this analysis, several event generator tunes are used. These are the PYTHIA6 (version 6.426 Ifl2l ) 
tunes Z2, Z2*, and CUETP6S1 HZl, PYTHIA8 (version 8.175 OH) tunes 4C OH and CUETP8S1 
IflTl . and HERWIG++ 2.7 with tune UE-EE-5C [14,19 j|. All of these tunes use the CTEQ6L1 1201 
parton density function. The energy dependence of pxo in these tunes is parameterised as 
Pto(\/s) = Fxo 1 x (\/s/Eo ) £ r where p^ F , Eo, and e are tune parameters summarised in table 
[l] These parameters were obtained by tuning to different data sets. The 4C tune was derived 
from early LHC data on charged particle multiplicities IfTHI . Z2, Z2*, CUETP6S1 and CUETP8S1 
were tuned to the previous UE results from CMS at 0.9 and 7 TeV |5l. In addition, CUETP6S1 
tune also included CDF data Ii2ll at 0.3, 0.9, and 1.96 TeV, while CUETP8S1 tune used 0.9 and 
1.96 TeV data for tuning. The UE-EE-5C was tuned with the ATLAS UE data at 0.9 and 7 TeV (7J 
and CDF UE data at 0.3, 0.9, and 1.96 TeV |2T| . None of the tunes make use of data at 2.76 TeV, 
making this a good test of interpolation between other centre-of-mass energies. The detector 
response was simulated in detail by using the GEANT4 package |22| , and simulated events were 
processed and reconstructed in the same manner as collision data. 
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Table 1: Summary of the parameters of the Monte Carlo generator tunes. 


Tune 

pV (GeV) 

Eo (GeV) 

e 

Z2 

1.832 

1800 

0.275 

Z2* 

1.921 

1800 

0.227 

CUETP6S1 

1.9096 

1800 

0.2479 

4C 

2.085 

1800 

0.19 

CUETP8S1 

2.1006 

1800 

0.2106 

UE-EE-5C 

3.91 

7000 

0.33 


3 Experimental methods 

3.1 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter, providing a magnetic field of 3.8 T. Within the solenoid volume there are several com¬ 
plementary detectors: a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic 
calorimeter, and a brass and scintillator hadron calorimeter, each composed of a barrel and 
two endcap sections. The silicon tracker measures charged particles within the pseudorapidity 
range |;/| < 2.5. For non-isolated particles of 1 < pj < lOGeV and ;/| < 1.4, the track reso¬ 
lutions are typically 1.5% in pj, and 25-90 (45-150) pm in the transverse (longitudinal) impact 
parameter |23] . Two of the CMS subdetectors acting as LHC beam monitors, the Beam Scintil¬ 
lation Counters (BSC) and the Beam Pick-up Timing for the experiments (BPTX) devices, are 
used to trigger the detector readout. The BSC are located along the beam line on each side of 
the Interaction Point (IP) at a distance of 10.86 m and cover the range 3.23 < |f/| < 4.65. The 
two BPTX devices, located inside the beam pipe at distances of 175 m from the IP, are designed 
to provide precise information on the bunch structure and timing of the incoming beams, with 
a time resolution better than 0.2 ns. A more detailed description of the CMS detector, together 
with a definition of the coordinate system used and the relevant kinematic variables, can be 
found in Ref. ffi24l . 

3.2 Event and track selection 

The present analysis is performed with a data sample of proton-proton collisions collected 
with the CMS detector at \/s = 2.76 TeV during a dedicated run in March 2011, corresponding 
to an integrated luminosity of 0.3 nb In 6.2% of the events there is an extra (pileup) pp 
collision, corresponding to an average of 0.12 overlapping pp collisions. Minimum bias events 
were recorded by requiring activity in both BSC counters in coincidence with signals from both 
BPTX devices (in contrast to Ref. ||fQ, where only one of the BPTX devices is required). To 
reduce the statistical uncertainty for the highly prescaled minimum bias trigger at large Pj , 
single-jet triggers based on information from the calorimeters, with pj thresholds at 20 and 
40 GeV, were also used to collect data (differently from Ref. |5|, where thresholds of 30 and 
50GeV are used). Events identified as originating from beam-halo background were removed 
from the sample l25l . The event selection requires exactly one primary vertex with more than 
four degrees of freedom (approximately 4 particles) and located no more than ±10 cm from the 
centre of the luminous region (beamspot) in the z-direction. 

For each selected event, the reconstructed track collection needs to be freed from undesired 
tracks, namely secondaries and background from track combinatorics and beam halo tracks. 
Tracks not corresponding to actual charged particles (misreconstructed tracks) are suppressed 
by imposing the high-purity selection criteria 123| . Secondary decays are suppressed by requir¬ 
ing that the impact parameter significance do/a(do) (measure of the distance between the track 
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and the primary vertex in the xy-plane) and the significance in the z-direction d z /cr(d z ) to be 
each less than 3. In order to remove tracks with poor momentum measurement, we require the 
relative uncertainty in the momentum measurement cr(pj)/pj to be less than 5%. The aver¬ 
age reconstruction efficiency for the selected tracks is about 85% and drops to 75% for tracks 
with pj ~ 0.5 GeV and // ~ 2, while the track misreconstruction rate is about 2%, increasing 
to about 8% for tracks with pj ~ 0.5 GeV and \t] % 2. Track efficiencies are determined by 
matching the generated level and reconstructed level tracks. 

The event hard scale and reference direction, for the identification of the UE sensitive region, 
are defined using leading "track jets" Il26l or charged-particle jets. The use of track jets makes 
the transition of leading tracks to leading jets more continuous and extends the pj coverage 
to larger values. These jets are reconstructed from tracks with pj > 0.5 GeV and |f/| < 2.5 
using the Seedless Infrared-Safe Cone (SISCone) |[27l algorithm with distance parameter of 0.5. 
Although anti-fcx Il28l is now the preferred algorithm at the LHC, the SISCone algorithm is 
chosen in this analysis for compatibility with previous results |5|. Furthermore, a comparison 
of the UE activity obtained at generator level using SISCone and anti -kj algorithms has been 
performed, finding differences of only a few percent for pj < 20 GeV. From all reconstructed 
track jets with \rj\ <2 and pj > 1.0 GeV, the one with the largest Pj 1 is selected. Only events 
containing at least one track jet fulfilling these criteria are considered for this analysis. Jets are 
reconstructed with a matching efficiency of 80% at Pj 1 ~ 1 GeV and up to 95% for pj 1 > 20 GeV. 
Trigger conditions are chosen to keep the trigger efficiency as uniform as possible and close to 
100%. For the p^ ranges in [1,25), [25,50), and [50,100) GeV, we use the minimum-bias and the 
two single-jet samples, respectively, corresponding to about 11M, 50k, and 23k selected events. 

3.3 Observables 

In this analysis we follow the same methodology as in the previous studies of the UE activity 
in events with a leading charged-particle jet, carried out at i/s = 0.9 and 7TeV K9. Charged- 
particle jets and charged particles produced at central pseudorapidity (\p\ < 2) with p T > land 
0.5 GeV, respectively, are used to study the UE properties. The direction of the leading charged- 
particle jet in the event is used to select charged particles in the transverse region defined by 
60° < |A<J>| < 120°, where A cp is the relative azimuthal distance between a charged particle 
and the leading jet. The UE is measured in terms of particle and Hpj densities, as a function of 
the leading pj\ which is used as an estimate for the hard scale of the interaction. The particle 
density ((N c h) / [A//A(A<J>)]) and Zpx density ((I 2Pt) / [A7/A(A<J>)]) are computed, respectively, 
as the average number of primary charged particles, and the average of their scalar pj sum, 
each per unit of ?/ and of A cp. 

As suggested in Ref. l29l , the transverse region can be studied in detail by separating — in¬ 
dependently for the particle multiplicity and for the pj sum — the 60° < Acp < 120° and the 
—120° < Acp < —60° ranges, and identifying the regions with higher and lower activities, re¬ 
ferred to as transMAX and transMIN, respectively. The two regions should have roughly equal 
activities for most events since the dominant production channel, two-jet production, is topo¬ 
logically symmetrical. In a three-jet topology, the transMAX side will capture the activity from 
the third jet. The difference between the measured densities in the transMAX and transMIN 
regions is called the transDIF density. The resulting particle and Zpx densities are expected to 
be sensitive to different components of the UE activity. 

Since the transMAX region contains the third jet, while both transMAX and transMIN regions 
receive contributions from MPI and beam remnants, the transMIN activity is sensitive to MPI 
and beam remnants, and the transDIF activity is sensitive to harder initial- and final-state radi- 



3.4 Corrections and systematic uncertainties 


5 


ation. The present approach extends the methodology employed in Ref. 15) 


3.4 Corrections and systematic uncertainties 


The UE observables ((Nch) / [Ai]A(A<p)] and (EPt) / [A;/ A(A<p)}) described in Section 3.3 are 
reconstructed from selected tracks, with pj > 0.5 GeV and // < 2, in the region transverse 
to the leading track-jet. These measured observables are corrected for detector effects and se¬ 
lection efficiencies to reflect the primary charged-particle activity using a 2-dimensional, it¬ 
erative unfolding technique H30l based on response matrices that correlate the generated and 
reconstructed level observables. These matrices are constructed from the generator level and 


reconstructed level UE and pj observables for PYTHIA6 Z2 events; this procedure accounts for 
detector effects and inefficiencies. The unfolding matrices are applied to a PYTHIA8 4C sam¬ 
ple to estimate the systematic uncertainties related to the correction procedure. These vary 
between 0.2% and 4%, depending on the observable, region and pj. 


Several other sources of systematic uncertainties may affect the results. These include the im¬ 
plementation of the simulation of the track and vertex selection criteria, tracker alignment and 
material content, background contamination, trigger conditions, and pileup contributions. The 
uncertainty in the simulation of the track selection is evaluated by applying various sets of se¬ 
lection criteria and comparing their effects on the data and on the simulated events. The impact 
parameter significance ranges are varied by one unit around the nominal window resulting in 
an effect on the densities of 0.6-4%. Replacing the high-purity selection by the simpler require¬ 
ment A/i a y ers > 4 and A/pi X ei i ayers > 2 for the silicon strip and pixel detector layers, respectively, 
has an effect of up to 0.8%. Varying the fraction of misreconstructed tracks by 50% affects the 
densities by 0.4-0.6%. The description of inactive tracker material in the simulation is ade¬ 
quate within 5% |4|, and increasing the material densities by 5% in the simulation induces a 
change in the observables of 1%. The effects of tracker misalignment, precision in the IP posi¬ 
tion, and dead channels, evaluated by varying the detector conditions in the MC simulation, 
are each found to change the results by 0.1-0.3%. The effect of varying the trigger and vertex 
efficiencies within their uncertainties, as well as the effect of pileup contributions, all lead to a 
negligible effect. 


Systematic uncertainties are largely independent of one another, but they are correlated among 
data points in each experimental distribution. They are added in quadrature to the statis¬ 
tical uncertainties and are shown in all figures. Systematic uncertainties mostly dominate 
the statistical ones, which are often smaller than the data points. Table [2] shows a summary 
of the systematic uncertainties as a range from the minimum to maximum values as they 
vary across region, observable and pj . The transMAX and transMIN regions tend to have 
a larger total systematic uncertainty than the other regions and the (£ pj) / [A;/A(Ac/?)] observ¬ 
able tends to have a slightly larger total systematic uncertainty by about 0.2% compared to the 
{N c h) / [A/ / A ( Ac/> )]observable. The total systematic uncertainty is large at low pj 1 and decreases 
to a minimum at pj ~ 3 GeV and then rises again up to a plateau for pj > 20 GeV. 


4 Results 

In Fig. [T| the (a) particle and (b) Hpj densities, after unfolding, are shown in the transverse 
region, relative to the leading charged-particle jet, as a function of pj\ A steep rise of the 
underlying event activity in the transverse region is seen up to pj* ~ 8 GeV, followed by a 
"saturation" (plateau-like) region, with nearly constant multiplicity and small E pj density in- 
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Table 2: Summary of the systematic uncertainties (in percentage) due to various sources. The 
values range from the minimum to maximum from all regions, observables, and across differ¬ 
ent Pj l values. 


Source 

Systematic (%) 

Unfolding procedure 

0.2-4 

Impact parameter significance 

0.6-4 

Fraction of misreconstructed tracks 

0.4-0.6 

Track selection 

0.1-0.8 

Material density 

1 

Dead channels 

0.1 

Tracker alignment 

0.2-0.3 

Interaction point position 

0.2 

Total 

1.9-5.8 


crease. In Fig. [2j the (left panes) particle and (right panes) Zy T densities after unfolding are 
shown as a function of pj in the transverse region with maximum and minimum activities 
(transMAX and transMIN), respectively. In the transMIN region, the amount of UE activity is 
roughly half that in the transMAX region. The p ] j dependences observed in the two regions 
are also quite different. At high pj, the distributions show a slow rise in the transMAX region, 
while for transMIN the flattening of the UE activity as a function of p ] j l is more pronounced. 
The corresponding distributions in the difference between the transMAX and transMIN regions 
(transDIF) are presented in Fig.|^ The particle and Zpx densities both show a rise with p^, and 
the plateau-like region above pj ~ 8 GeV—seen for distributions in the individual transMAX 
and transMIN regions—is replaced by an increase as a function of p^. 



0 10 20 30 40 50 60 70 80 90 100 

pf [GeV] 



0 10 20 30 40 50 60 70 80 90 100 

pf [GeV] 


Figure 1: Measured (left) particle density, and (right) Zpx density, in the transverse region 
relative to the leading charged-particle jet in the event ( \rj\ < 2, 60° < Ac/; < 120°), as a 
function pj\ The data (symbols) are compared to various MC simulations (curves). The ratios 
of MC simulations to the measurements are shown in the bottom panels. The inner error bars 
correspond to the statistical uncertainties, and the outer error bars represent the statistical and 
systematic uncertainties added in quadrature. 


The rapid increase of the UE activity with Pj 1 in the region below ~8 GeV is mainly attributed 
to the increase of MPI activity as the hard scale of the interaction increases IfTTI . This fast rise 
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Figure 2: Measured (left panes) particle density, and (right panes) Zpx density, in the trans- 
MAX and transMIN regions (60° < |A<p| < 120°, relative to the leading charged-particle jet in 
the event, with maximum/minimum UE activity), as a function of The definitions of the 
symbols and error bars are the same as for Fig. [Tl 
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4 Results 
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CMS TransDIF density Vs = 2.76 TeV 



Figure 3: Measured transDIF activity (see text for its definition) for (left) particle density, and 
(right) Zpx density, as a function of Pj l . The definitions of the symbols and error bars are the 
same as for Fig. [lj 


is followed by a saturation region (for the transverse and especially transMIN distributions), 
with nearly constant multiplicity and small Zpx density increase. This behaviour is expected 
as a consequence of a nearly full overlap of the colliding protons in interactions yielding the 
hardest parton-parton scatterings. When pp collisions occur for very small impact parameter, 
the amount of MPI activity saturates [fTOl HTI . Such a distinct Pj ^dependent pattern in the 
amount of UE activity (sharp rise followed by a plateau above the pj* ~ 8 GeV transition) is 
clearly seen for all the observables presented, especially in the transMIN region. In contrast, 
the transMAX and transDIF distributions show a continuous rise with p ] j l also in the high-px 
regime. This is expected to be caused by contributions from initial- and final-state radiation 
in the transverse region Il29| . Following such an interpretation, the present results provide 
constraints on the modelling of the different UE components. 

The results are compared to recent tunes of the PYTHIA and HERWIG++ event generators. All 
PYTHIA6 and PYTHIA8 tunes predict the distinctive change in the amount of activity as a func¬ 
tion of the leading jet pj within 5-10%. The HERWIG++ UE-EE-5C tune also provides a fair 
description of the data. In general, the data-model agreement improves for the transDIF den¬ 
sities. The continuous increase observed at high-p^ 1 in the transDIF distributions is well re¬ 
produced by all MC tunes, corroborating the hypothesis of increased contributions of QCD 
radiation from the hardest scattered partons. The same trend is observed in pp collisions at 
1.96 TeV El. The latest PYTHIA6 (pythia8) tune CUETP6S1 (CUETP8S1) improves the descrip¬ 
tion of the data in comparison to the results obtained with the parameters of the previous Z2* 
(4C) tune. 

The centre-of-mass energy dependence of the UE activity in the transverse region is presented 
in Fig. [I] as a function of Pj 1 for -^/s = 0.9, 2.76, and 7 TeV 010- A fast rise with increasing 
centre-of-mass energy of the activity in the transverse region is observed for the same value 
of the leading charged-particle p^. This is expected from the higher parton densities probed 
at low-x in the protons, and the larger phase space available for parton radiation. All tunes 
predict a centre-of-mass energy dependence of the UE activity which is consistent with that of 
the data. 
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Figure 4: Comparison of UE activity at i/s = 0.9, 2.76, and 7TeV for (left) particle density, and 
(right) Zp T density, as a function of pj 1 fH|5l. The data (symbols) are compared to various MC 
simulations (curves). The definition of the error bars is the same as for Fig. [l] 

The measurements presented here provide constraints for the development and tuning of the 
underlying event description implemented in MC models. In particular, they may allow im¬ 
proving the modelling of key ingredients—such as multiparton interactions, QCD radiation, 
energy evolution of the transverse proton profile, etc.—, which will play an increasing role at 
higher proton-proton collision energies. 


5 Summary 

The measurement of the underlying event (UE) activity in proton-proton collisions at y/s = 
2.76 TeV has been presented using events with a charged-particle jet produced at central pseu¬ 
dorapidity (| pi et | < 2) with transverse momenta 1 < Pj < 100 GeV. This analysis comple¬ 
ments the results of previous similar measurements at y/s = 0.9 and 7 TeV. 

The UE activity is measured in the transverse region and further studied in terms of the trans- 
MAX, transMIN and transDIF activities. A steep rise of the underlying activity in the trans¬ 
verse region is seen with increasing leading jet pj. This fast rise is followed by a leveling above 
Pj* m 8 GeV, with nearly constant particle density and small E pj density increase. Such a dis¬ 
tinct pattern (fast rise followed by a leveling of the UE hadronic activity) is clearly seen for all 
the observables in the various regions, and is compatible with the impact parameter picture of 
pp collisions featuring an increasing number of MPI for increasing overlap followed by a satu¬ 
ration of hadron production once the hardest most-central collisions are reached. The transDIF 
density distributions show an increase of the activity as a function of pj, corroborating the 
hypothesis of more intense ISR and FSR from the increasingly harder parton-parton scatter. 

The results are compared to recent tunes of pythia and HERWIG++ Monte Carlo event gener¬ 
ators. The PYTHIA6, PYTHIA8, and HERWIG++ tunes describe the data within 5 to 10%. All MC 
tunes predict a collision energy dependence of the hadronic activity similar to that observed in 
the data. The ability of the latest Monte Carlo generator tunes to describe the data confirms the 
validity of the tunes and lends confidence to the predictions of UE activity for higher collision 
energies. 
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